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TECHNICAL NOTE 3839 

EXPERIMENTAL DROPLET IMPINGEMENT ON SEVERAL TWO-DIMENSIONAL 
AIRFOILS WITH THICKNESS RATIOS OF 6 TO 16 PERCENT 
By Thomas F. Gelder, William H. Smyers, Jr., and Uwe von Glahn 


SUMMARY 

The rate and area of cloud droplet impingement on several two- 
dimensional swept and unswept airfoils were obtained experimentally in 
the NACA Lewis icing tunnel with a dye-tracer technique. Airfoil thick- 
ness ratios of 6 to 16 percent, angles of attack from 0° to 12°, and 
chord sizes from 13 to 96 inches were included in the study. The data 
were obtained at 152 knots and are extended to other conditions by dimen 
sionless impingement parameters. 

In general, the data show that the total and local collection effi- 
ciencies and impingement limits are primary functions of the modified 
inertia parameter (in which airspeed, droplet size, and body size are 
the most significant variables) and the airfoil thickness ratio. Local 
collection efficiencies and impingement limits also depend on angle of 
attack. Secondary factors affecting impingement characteristics are air 
foil shape, camber, and sweep angle. The impingement characteristics ob 
tained experimentally for the airfoils were within ilO percent on the 
average of the characteristics calculated from theoretical trajectories. 
Over the range of conditions studied, the experimental data demonstrate 
that a specific method can be used to predict the impingement character- 
istics of swept airfoils with large aspect ratios from the data for un- 
swept airfoils of the same series. 


INTRODUCTION 

Knowledge of the local and total rates of cloud droplet impingement 
and of the surfacewise extent or limit of droplet impingement on bodies 
is required for the design and evaluation of icing -protect ion equipment 
for aircraft. These impingement characteristics are important factors 
in determining the extent of the surface to be protected, the shape and 
location of some ice formations on aircraft components, the aerodynamic 
penalties associated with icing of aircraft surfaces, and the local and 
total requirements for various thermal and fluid protection systems. 
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Previous studies (refs. 1 to 12) report the droplet trajectories 
about several two-dimensional bodies and bodies of revolution. These 
studies used differential analyzers for computing the droplet paths after 
the flow field about the body had been obtained. An empirical method for 
obtaining the impingement characteristics of airfoil sections is presented 
in reference 13. This method, however, is more suited to airfoils with 
blunt leading edges, because the basic data used in developing the method 
were obtained from four Joukowski airfoils and only one low-drag airfoil. 
For the two-dimensional case, a method for applying trajectory data from 
unswept airfoils to swept airfoils is presented in reference 14. 

Droplet trajectories about bodies with unknown or complex flow fields 
are difficult to obtain with a differential analyzer. Therefore, a wind- 
tunnel method using a dye -tracer technique to obtain experimentally the 
impingement characteristics of bodies has been developed (ref. 15). In 
this technique water treated with known small quantities of a water- 
soluble dye is sprayed into the tunnel airstream by nozzles a large dis- 
tance ahead of the body. The surface of the body is covered with blotter 
paper or a similar absorbent material upon which the dyed droplets im- 
pinge and are absorbed essentially upon contact. At the point of droplet 
impact and absorption, a permanent dye trace is deposited. The amount of 
dye deposited in a measured time interval can be determined by a colori- 
metric analysis of the blotter paper and can be converted into the amount 
of impinged, water that produced the dye trace. From such an analysis and 
from known values of spray-cloud water content and droplet sizes, the im- 
pingement characteristics of a body can be determined readily, as dis- 
cussed in reference 15. 

In an extensive program of icing studies conducted in the NACA Lewis 
icing tunnel on various two- and three-dimensional bodies, experimental 
impingement data on six swept and two unswept airfoils (all two- 
dimensional) have been obtained. Although the airfoils used in these 
studies were a rather ad hoc collection of shapes and sizes, this report 
makes these data generally available and correlates the data as much as 
possible. The impingement data were obtained with airfoil chord sizes 
ranging from 13 to 96 inches, three volume-median droplet sizes ranging 
from 11 to 19 microns in diameter, and a nominal airspeed of 152 knots. 

The thickness ratio of the airfoils studied varied from 6 to 16- percent. 

The airfoil impingement rates and limits obtained are presented in 
terms of dimensionless impingement parameters. These dimensionless pa- 
rameters allow interpolation and extension of the experimental results 
over a wide range of operating conditions. The experimental impingement 
values for several airfoils are compared with values calculated from 
theoretically obtained trajectory data. 
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APPARATUS 
Airfoil Models 

This study of droplet impingement on various airfoil sections was 
conducted in the 6- by 9 -foot test section of the MCA Lewis icing tun- 
nel. The models, unless otherwise noted, were made of wood and spanned 
the 6-foot height of the tunnel (fig. l). The airfoils are listed in 
the following table, and dimensionless streamwise sections are presented 
in figure 2. 


Airfoil section 
(fig. 2) 

Chord length 
in streamwise 
direction, in. 

Remarks 

(a) Joukowski 0015 

13 


(b) Joukowski 0015 

96 

Smooth sheet -metal surface 

(c) 632-015 

13 


(d) 652-015 

13 


(e) 65,2-216 

96 

Isheet -metal surface modified by 
3/l6-in. -thick de-icing boot extend- 
ing from s^ of 0.156 to s^ of 

0.250 (ref. 16) 

(f) 65j_-212 

13 


(g) 65^-212 

72 


(h) 65 ;l“206 

72 

Although the low-drag range for this 
series airfoil is <0.1 for thickness 
ratios <0.12 and thus the subscript 
1 is usually omitted, it is retained 
herein to preserve similarity with 
the 65]_-212 section 

(i) 65i-212 

87.9 

Swept 35^^, design section in plane 

(j) 65 ;l“206 

87.9 

perpendicular to leading edge 


^Symbols are defined in appendix A. 


The leading edge of the models was about 1.5 and 2.2 chord lengths from 
the entrance of the test section for the 96- and 72-inch-chord airfoils, 
respectively, and 9 chord lengths for the 13-inch-chord airfoils. These 
longitudinal locations indicate the length of the upstream flow field. 
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The local pressures on the surface of several models were obtained 
"by use of pressure belts. These experimental data, uncorrected for 
tunnel wall effects, were used to indicate variations from the theoret- 
ical pressure distributions. 


Spray System and Related Equipment 

The spray cloud was provided by air-water atomizing nozzles located 
in the quieting chamber upstream of the tunnel test section. The noz- 
zles were always positioned to provide a cloud that was relatively uni- 
form in liquid-water content and droplet-size distribution in the test 
section. The dye-water solution and air pressures to the nozzles were 
set by means of pressure transmitters and manometers. The spray was 
turned on and off by fast -act ion solenoid valves, while the spray dura- 
tion was set and recorded by an electric timer. Further details of the 
spray system are described in reference 15. 


PROCEDURE 
Blotter Mounting 

For the larger models, a 3-inch-wide blotter was rubber- cemented to 
a vellum strip, which in turn was cemented to the airfoil surface as 
shown in figure 1. The cementing prevented the blotter from being lifted 
from the airfoil surface by aerodynamic forces. The edges of the blotter 
were also taped to the airfoil surface. The vellum strip prevented dam- 
age to the blotter during removal from the model. After exposure to the 
spray cloud the vellum and blotter were removed as a unit and later 
separated carefully. For the small models a 2-inch-wide blotter was 
stretched tightly over the leading edge and taped to the airfoil surface 
along all the blotter edges. 


Tunnel Conditioning and Blotter Exposure 

In order to minimize the evaporation of the droplets during their 
time of travel from the spray nozzles to the body (about 3/4 sec), the 
entire tunnel airstream was nearly saturated before the body was ejqposed 
to the dyed-water spray. Saturation of the test-section airstream was 
achieved through the control of tunnel air temperature and the addition 
of steam into the tunnel until a light condensation cloud resulted. 

The studies reported herein were conducted at the following nominal 
conditions: Free-stream velocity, 152 knots; static pressure, 28.1 

inches of mercury; and static air temperature, 50° F. ■ 
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The procedure for each run was to preload the air and dyed-water 
pressure in the spray system and preset the exposure time. (Air-water 
gage pressure ratios similar to those used in ref. 15 - i.e.^ 0.5^ 0.6, 
and 0.8 - were used herein to obtain impingement data; low pressure ratio 
was used to obtain large droplet sizes, while high pressure ratio was 
used to obtain small droplets.) With the tunnel air properly conditioned 
as to speed, temperature, and humidity, the blotter -wrapped model was ex- 
posed to the dyed spray for the preset time interval. After tunnel shut- 
down, the blotter was removed from the model. 

In these studies the exposure time varied from 2 to 7 seconds for 
the 13-inch-chord airfoils and from 3 to 12 seconds for the large -chord 
airfoils at air-water pressure ratios of 0.5 to 0.8, respectively. 

For each air-water pressure ratio, a relatively uniform cloud with 
local liquid-water-content variations within £10 percent and essentially 
the same droplet-size distribution were obtained. The reproducibility 
of the average liquid-water content from one model exposure to the next 
was about ±5 percent. 


Spray Cloud Properties 

The cloud total liquid-water content was obtained by collecting dye 
from the spray cloud in an aspirating device (a tube that draws in air 
and liquid water at free-stream conditions, ref. 15) . The inlet velocity 
of the device was always within 1 percent of the free-stream value, denot- 
ing theoretically a 100-percent collection efficiency. 

The droplet -size distribution was determined by the method outlined 
in reference 15, in which the experimental impingement rates for cylinders 
are related to theoretical data for similar cylinders. In the present 
study, however, small 36. 5 -per cent -thick Joukowski airfoils were used 
instead of cylinders. The absolute values of droplet size from this air- 
foil section (see appendix B) generally confirm those obtained from the 
cylinders of reference 15, but the body size trend is reduced from that 
of reference 15. 

The ratio of droplet diameter to voliame -median diameter^ is presented 
in figure 3 as a function of the ratio of cumulative liquid-water content 
to total liquid-water content. The volume-median droplet diameters are 
believed accurate within £6 percent. These data are from the aspirator 
and 36. 5-percent Joukowski airfoil analyses for the three spray conditions 
(air-water pressure ratios) studied. 


Volume-median droplet diameter is that diameter for which half the 
total liquid-water content is contained in droplets larger than the vol- 
ume median and half in droplets smaller than the volume median. 
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The pertinent spray cloud properties are suimarized in the follow- 
ing table; 


Air- 

water 

pressure 

ratio 

Air 

pressure, 
Ib/sq in. 
gage 

Water 
(dye so- 
lution) 
pressure, 
Ib/sq in. 
gage 

Approx, 
max. 
droplet 
diam. , 

^ax^ 

microns 

Range of 

total 

liquid- 

water 

content, 

^t? 

g/cu m 

Volume- 
median 
droplet 
diam. , 

^ed^ 

microns 

0.5 

60 

120 

59 

0.46-0.65 

18.6 

.6 

60 

100 

48 

.37- .50 

16.7 

.8 

80 

100 

29 

.22- .33 

11.5 


The number and spacing of the spray nozzles varied during the course of 
the airfoil program because of other unrelated test programs interspersed 
between those reported herein. These changes resulted in a range of 
cloud liquid-water content, as noted in the previous table, and were ac- 
counted for in analyzing the experimental impingement data. Droplet-size 
distribution was not affected by these nozzle changes. 


Colorimetric Analysis 

In the colorimetric analysis of the dyed blotter, small 
^ ^ ^ 1-in.^ segments of area AAg are punched from the 

blotter as shown in figure 4. The dye is dissolved out of each segment 
"with a known quantity of distilled water (ref. 15). The concentration 
of this solution is determined by the amount of light of a suitable wave- 
length transmitted through the solution in a calibrated colorimeter. The 
amount of dye collected on the segment is converted into the weight of 
water (dye-water solution) that impinged on_the blotter segment during 
the exposure. The local impingement rate for a segment as given in 

reference 15 is expressed as 


77 0.794 Fb 

“ t^g'AAs ' 


lb water 



( 1 ) 


ANALYSIS OF DATA 

The analysis of the data obtained from the dye -impingement records 
consists in evaluating the local and total collection efficiencies of 
the airfoils and the extent or limit of impingement on the airfoil sur- 
faces. In order to analyze the experimental data, the water content and 
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droplet-size distribution of* the spray cloud also must be known. Final- 
ly, methods of extending the data for conditions other than those used 
in the tests must be employed in order to render the data generally use- 
ful. A detailed discussion of the dye analysis is presented in reference 
15 and reviewed herein for convenience in presenting the experimental 
data. 


The local rate of water-droplet impingement and limit of im- 
pingement are obtained as a direct result of the dye-tracer tech- 

nique used herein. The dimensionless impingement parameters ^ and 
are obtained from the following equations (ref. 15): 


P = 


W, 


0.329 UqW^ 


( 2 ) 


_ J 

E. 


^s 


u,max 


Wp dAc 


"m 


’Zj_max 


0.329 UqW^A-p A-pO Si TTiay 


i , m 


max 


P dA, 



For a two-dimensional airfoil, equation (3) is rewritten for con- 
venience as 


u,max 


Em = 


Wp ds 

Z,max 1 I 

1.329 UQW-f-h ” h 'J st „ 


max: 


P ds 


Z,max 


(4) 


Equation (4) is based on projected frontal height h rather than on the 
airfoil thickness ratio used by some investigators. Figure 5 shows ratio 
of- projected frontal height to chord length h plotted against angle of 
attack a for the airfoils used herein. For the swept airfoils, h and 
s are referred to the free-stream direction. 


Total collection efficiency and impingement limits are often pre- 
sented in terms of K and cp, where K indicates the inertia of the 
droplet and <p represents the deviation of the droplet drag forces from 
Stokes' law, for correlating impingement characteristics. Reference 17 
discusses and illustrates previously determined analytical airfoil- im- 
pingement data in terms of a modified K parameter defined as 
Kq = K(X/Xg). The ratio X/Xg is a function of Rcq as shown in fig- 
ure 6 (data from ref. 12). Plotting a dependent impingement parameter 
such as Eqj or s ^^^y . as a function of Kq yields a set of experiment- 
ally or analytically determined points that can be essentially represented 
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by a single curve, independent of cp . This curve is approximately the 
solution obtained using Stokes’ law for sphere (droplet) drag. The ex-> 
tention or interpolation of experimental as well as analytical data 
points over a wide range of the pertinent impingement variables (droplet 
size, body size, speed, and altitude) is greatly facilitated by this Kq 

parameter, even though no complete theoretical proof of its significance 
or validity is available at this time. 

The impingement parameters previously discussed are often presented 
in the literature in terms of clouds containing droplets all of the same 
size. Analytical calculations show that the presence of a droplet-size 
distribution does not alter the usefulness of the Kq parameter if it is 

evaluated with the following droplet sizes: (l) The use of volume-median 
droplet size to calculate a "weighted" Kq will reduce weighted total 

collection efficiency to data representable by a single curve; and (2) the 
use of maximum droplet size in calculating Kq will reduce limits of 

impingement to data representable by a s_ingle curve. Reasonable exten- 
sion and interpolation of experimental Eni and Smax data obtained with 
droplet-size distribution for conditions other than those studied are 
possible, therefore, with the or Kq parameter, respectively. 

Correlation of 3 (uniform droplet size) or p* (distribution of drop- 
let sizes) with was possible only if each value of p or ^ 

used was obtained for the same numerical value of s’, where s’ is the 
dimensionless surface distance_measured from the s location of p or ^ 
to the location of P^iax ^max.» I’espectively. 

From the theoretical and experimental impingement results, it was 
determined that the surface location of Pj^gx ^max^ measured from the 

zero-chord point, denoted as s", does not occur at the same s location 
for various Ko,med 'values except for symmetrical airfoils at zero angle 

of attack. Generally, s" occurs between the air stagnation point on the 
airfoil (max. pressure point) and the foremost point on the airfoil. The 
foremost point on the airfoil is where the airfoil surface is perpendicular 
to the free-stream-velocity direction. As the value of the modified iner- 
tia parameter increases, s" moves toward the foremost point on 

the airfoil, because the droplet paths approach straight-line trajectories. 
As the value of decreases, s" moves toward the maximum pressure 

point, because the droplet inertia is approaching that of air particles. 

The following empirical method of analysis was adopted for correlat- 
ing p with KQ^^e, 3 ^. From plots of the experimental data of p" against 
s, values of p were selected at specified values of s’. Plotting F 
as a function of KQ^med various s’ values yields points reasonably 
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represented by a single curve. The relation of these (3 values to their 
true surface location s is accomplished by a plot of s” as a function 
of KQ^med relation 


s = s’ + s" ( 5 ) 

A negative sign herein denotes the airfoil upper-surface values; a 
positive sign denotes the airfoil lower surface. 


RESULTS 

A complete tabulation of the local collection efficiency p”, the 
maximum extent of impingement on the airfoil surface s^^^, and the total 

collection efficiency for each airfoil and impingement condition 

studied is listed in table li The p values are tabulated as a func- 
tion of s, the surface distance from the zero -chord point on the airfoil 
divided by chord length. These p data_are calculated from faired val- 
ues of Wp and equation (2). Typical values as a function of s 

are shown in figure 7 for several repeat runs. The data in this figure 
show a repeatability of better than ±10 percent. 

In order to emphasize the significant trends and variables affecting 
the impingement characteristics of airfoils, this section of the report 
presents the experimental data in terms of (l) typical p curves as a 
function of s and (2) dimensionless Kg^med ^0 max parameters. 

Part (l) consists in a. general evaluation of the effect on impingement 
characteristics of the basic airfoil geometry, including such items as 
airfoil angle of attack, camber, thickness ratio, airfoil shape, and 
sweep angle. Part (2) presents the effects on airfoil impingement char- 
acteristics of varying the droplet size, airspeed, and model size as ex- 
pressed by a variation of the dimensionless Kq and Kq 

parameters. 


Effect of Airfoil Geometry on Impingement Characteristics 

In order to illustrate the effect of airfoil geometry and attitude 
on typical values of local collection efficiency, p is presented as a 
function of surface distance divided by chord length s (fig. 8) . The 
data used in figure 8 were obtained from table I. 

Angle of attack . - An increase in angle of attack for the same spray 
cloud conditions will increase the extent of impingement on the lower 
surface of an airfoil decrease the extent on the upper sur- 

face s■^^J^a,x• figure 8(a) (Joukowski 0015 airfoil at angles of attack 
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of 4^, and 8*^) the impingement limit on the lower surface increased 
from a value of 0.158 to 0.379 as the angle of attack was 

increased from 0° to 8°. Concurrently, on the upper surface s^ 
decreased from -0.158 to -0.058. 

Generally, for a given s the local p values on the lower surface 
increase while those on the upper surface decrease as the angle of attack 
is increased. The value of ^ tendency (both from experimental 

and theoretical data) to decrease slightly (less than 10%) as the angle 
of attack is increased from 0° to 8°. As the angle of attack is in- 
creased, the v-alue is located farther aft along the lower surface, 

as shown in figure 8(a). At 8° angle of attack occurs at s" of 

0.012 as compared with £" of 0.005 and 0 at 4*^ and 0°, respectively. 

The shape of the local p curve is symmetrical or nearly symmetrical 
(depending on whether the airfoil is symmetrical or is cambered, respec- 
tively) at zero angle of attack. As the angle of attack is increased, 
the p curve has a steeper impingement gradient on the upper surface and 
a lesser gradient on the lower surface. For the airfoils, droplet sizes, 
and operating conditions used herein, the location of Py.,0.^ occurs be- 

tween the air stagnation region and the foremost point of the airfoil. 

The values for the 12- to 16 -percent thick airfoils show no 

great change with an increase in angle of attack (table I) over the range 
of values studied (0.0057 to 0.093). However, these airfoils 

with increasing angle of attack will have an increased total water catch 
per foot of span almost proportional to the increased projected frontal 
height h of the airfoil (as will be discussed later). The total collec- 
tion efficiency for the NACA 65p-206 airfoil also does not vary apprecia- 
bly for the limited range of angles of attack and Kq .values studied 

(0° to 4.3° and 0.0077 to 0.0167, respectively). According to theoretical 
data for a thin airfoil (ref. 11), however, the Em for the NACA 65A-004 
airfoil increased markedly between angles of attack of 0*^ and 8°, the Em 
being 40 to 80 percent greater at 8° than at 0° for Kq values of 0.01 
and 0.10, respectively. Therefore, the Em for the NACA 65i“206 may in- 
crease for angles greater than 4° in the range of 0.01<Kq ^^^^<0.10 in a 
manner similar to that for the NACA 65A-004. 

I 

Airfoil thickness . - A change in airfoil thickness ratio (for the 
same airfoil series) has mixed effects. An increase in thickness ratio 
causes a decrease in as shown in figure 8(b) (NACA 65]_-206 and 

65i-212 airfoils both at 0° angle of attack). The value of for 

the thin 65 q_- 206 airfoil is 27 percent greater than that for the thicker 
651-212 airfoil. Between the region of Pmax limits of impinge- 

ment, local p values for a thick airfoil are generally higher than those 
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of a thin one. In figure 8(h), the upper-surface limit of impingement 
Su,max on the 12 -percent -thick airfoil is -0.06, while that on the 6- 
percent-thick airfoil is -0.09. The lower-surface limit of impingement 
sz,max is 0.03 for the thick airfoil compared with 0.02 for the thin 
airfoil. For the range of cloud properties, angle of attack, and airfoil 
series studied herein, the higher values of ^^lgQc smaller h of 

the 651-206 airfoil yielded a higher total collection efficiency than 
those of the 65i-212 airfoil. (See table I, eq,. (4), and subsequent 
discussion. ) 

Airfoil shape . - The MCA 652-015 and 632-015 and the Joukowski 0015 

airfoils are compared in the same cloud conditions and at 0*^ angle of 
attack in figure 8(c). These three airfoils are synmetrical and have a 
maximum thickness of 15 percent; they differ in the location of maximum 
thickness and leading-edge radius (see fig. 2). The Joukowski 0015 is 
the bluntest, with the maximum-thickness point at 25-percent chord; the 
652-015 is the sharpest, with the maximum-thickness point at 40-percent 

chord; and the 632-015 airfoil is intermediate, with maximum thickness 
at 35-percent chord. The blunt Joukowski 0015 has a lower P^ax ^silue 

but higher local P values farther aft on the surface than the sharp, 
low-drag 652-015 airfoil. The value of 652-015 airfoil 

was about 20 percent greater than for the Joukowski 0015. The 632-015 
airfoil data show P^ax values between those obtained for the other two 

15 -percent -thick airfoils. In addition, the limit of impingement on the 
bluff airfoil is less than that on the sharp low-drag airfoil, as shown 
in figure 8(c) . 

The total collection efficiency of the blunt airfoil (Joukowski 0015) 
in the range of KQ^med covered herein (0.04 to O.l) is 20 to 40 percent 

higher than that of the low -drag airfoil (652-015) of the same thickness 
ratio (table l). In the determination of the values for these air- 
foils, the higher F^ax greater total impingement area of the low- 

^ag airfoil (fig. 8(c)) are insufficient to offset the generally higher 
p values (except near P*niax) blunt airfoils. 

Airfoil camber . - The extent of impingement at 0° angle of attack 
may be greater on the upper surface than on the lower surface for a cam- 
bered airfoil (fig. 8(b)); however, this is not generally true for angles 
greater than 0°. For the MCA 65;^-212 airfoil, the maximum impingement 

limit on the upper surface s^^max is -0.06, whereas on the lower surface 
the impingement limit is 0.03. No direct comparison of the effect of air- 
foil camber on local or total collection efficiency can be made with the 
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data available_, because the same airfoil series vas not studied vith and 
■without camber. An interpolation, however, of the Em data for the 
cambered 65-series airfoils (653_-206, 65]_-212, and 65,2-216, the last 

airfoil having similar geometric shape to the other two, 1)) to a 

cambered 15 -percent -th^ck airfoil can be made. Comparison of these in- 
terpolated values of E^ with those of the uncambered CSg-OlS airfoil 

indicates that the small camber involved has only a secondary effect on 
the total collection efficiency. 

Airfoil sweepback . - Sweeping back an airfoil and keeping the same 
physical shape (yawing the airfoil) generally have only a small effect 
on impingement limits and local and total collection efficiencies . Com- 
paring the data for the unswept and swept 65-series airfoils (table I) 
shows that the collection efficiencies and limits of impingement are a 
little less, in general, on the yawed or swept airfoils than on the un- 
swept airfoils. A correlation of swept- and unswept -airfoil impingement 
data is presented later in the DISCUSSION. 


Dimensionless Presentation of Data 

The experimental impingement characteristics of airfoils E^, 
and are conveniently presented as a function of the pertinent 

modified inertia parameter or Kq for purposes of extra- 

polation and comparison in figures 9 to 12. Presentation of data in 
this form permits a ready evaluation of airfoil impingement characteris- 
tics in terms of droplet size, air temperature, altitude, component size, 
and airspeed. 

Variation of p with ” In. figure 9, average p values at 

selected surface locations s’ are shown as functions of Kq for 

all airfoils studied. These locations of 3 are referenced to the loca- 
tion of P|nax discussed previously and were obtained from the origi- 
nal curves of 3 against s. It is apparent from the curves shewn in 
figure 9 that the method of analysis for F (outlined in the ANALYSIS) 
will produce data that can be essentially represented by a single curve 
for particular s’ locations. 

In the range of 0.04 < <0.1 (l3-in. -chord airfoil data) the 
point value of 3 ma x could not be readily ascertained because of the 
width of the blotter punch used. Generally, the 3max curve for the 

13-inch-chord airfoils is estimated from the more reliable data in the 
range of 0.001 < <0.01 and from the shape of the 3 curves aft of 

^max ~ 0.01, etc.) for the values of greater than 0.04. 
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The local collection efficiency p increases with an increase in 
^0,med shown in figure 9. The order of magnitude of the increase in 

^ ^0,med increases depends on the surface location and is a complex 

function of the airfoil ^gle of attack and shape. In the range of 
0.01< <0.1, the p values are believed accurate to ±10 percent; 

but, for Pmax ^,med Sweater than 0.02, a ±25 -percent accuracy is 

estimated. 

Ui 
CO 
CO 

to 


Variation of s^^g^ with KQ^^ia^ . - An increase in Kg^max ^ili in- 
crease the extent of impingement on both surfaces of the airfoil as shown 
in figure 11. With increasing there Is a greater change in the 

lower-surface limit of impingement on a low-drag airfoil (MCA 652_-212) 

than on a blunt airfoil (joukowski 0015). On the upper surface no maxked 
trends with airfoil shape are apparent except at 0° angle of attack. As 
previously discussed (fig. 8(a)), an increase in angle of attack will in- 
crease the extent of impingement on the lower surface and decrease the 
extent on the upper surface for all airfoils over the entire range of 
^0,max v-alues. 

Variation of with - The total collection efficiency 

of an airfoil increases as increases (fig. 12). These data 

are plotted from table I for gm angle of attack of 0°. In the range of 
^0,med'^^*^^ increase in Ej^ with a decrease in the thickness ratio 

is readily apparent in figure 12 for airfoils of the same series and 
camber. For example, at 0° angle of attack and Kg of 0.01 the 

values for the NACA 65i-206, 65^-212, and 65,2-216 airfoils are 0.105, 

0.06, and 0.05, respectively. These airfoils, although their numbering 
systems are somewhat different, are of the same series, differing prima- 
rily in thickness ratio and only to a minor degree in leading-edge radius 
and location of maximum-thickness point. These latter differences are 
considered of secondary significance in the evaluation of total collection 
efficiency. Reversals of this change in E^j with thickness ratio may 
occur at high values 0.02 ), especially as the thickness 

ratio approaches 6 or 4 percent. 


The curves of figure 10 present s” (the surface location of 6 ) 

^TriB,X ^ 

as a function of Because of the width of the blotter punch 

used, s’* could not be precisely established. Consequently, s" is 
represented in figure 10 by a dashed line. 
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The preceding E-^ comparisons do not necessarily mean a similar 

comparison for total water caught on an airfoil. For a two-dimensional 

wing the total water catch per foot span is 

= 0.329 UQW^ch (6) 

Total catch is thus proportional to the product of projected frontal 

height h, and chord length c. For example^ at 0° angle of attack and 
Kq of 0.007, the NACA 65,2-216 airfoil has an E^^ value of 0.032 

and a value proportional to 0.16x0.032 = 0.0051 (values from figs. 

5 and 12). In comparison, the NACA 65]_-206 has a higher Ejj^ value of 
0.085 hut Wm proportional to 0.06x0.085 = 0.0051 is the same. A similar 
comparison at Kq of 0.013 makes proportional to 0.0109 and 

•0.0072 for the NACA 65,2-216 and 652^-206, respectively, although E^^ for 

the thicker airfoil is about 40 percent less than that for the thinner 
one . 

The effect of angle of attack on E^ for all the airfoils studied 

at one spray condition is presented in figure 13. These data were ob- 
tained at a free-stream velocity Uq of 152 knots and a volume-median 

droplet diameter d^^^^ of 16.7 microns (air-water pressure ratio, 0.6) 

and are typical of the other spray cloud conditions.^ (Cross plots of 
these data, obtained from table I, yield curves of Ej^ against 

similar to those of fig. 12.) As previously discussed, Em for the air- 
foils and spray conditions studied does not vary appreciably with angle 
of attack. 

Because the h difference among the airfoils studied lessen as 
angle of attack increases (fig. 5), values of total water catch per foot 
span (see ep. (6)) will depend largely on the values of Ejjj^ for 

angles of attack greater than 4^. 

DISCUSSION 

The following discussion is based on comparison of the experimental 
results with available theoretical trajectory data. In order to compare 
experimental impingement data with that calculated from trajectory data 
for the same body, the flow fields in both cases must be similar. The 
local velocity distributions obtained experimentally on several of the 
airfoils are shown in figure 14. Also shown in this figure are the theo- 
retical values used to set up the flow field for the trajectory calcula- 
tions of reference 2. In general, the experimental results agree well 
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(especially near. the leading edge) with the theoretical. The experi- 
mental velocity distributions are generally slightly high on both sur- 
faces by an average of 5 percent. At angles of attack up to 4°, the 
local velocity data from the large-chord airfoils agrees well with those 
for the small-chord airfoils. 


Impingement Characteristics 

The theoretical impingement characteristics for the Joukowski 0015 
airfoil at angles of attack of 0° or 4° and the unswept NACA 65j_-212 air- 
foil at an angle of attack of 4° can be obtained from the trajectory 
studies of references 5 and 2, respectively. These trajectory results, 
when "weighted" (as described in ref. l) for the droplet-size distribu- 
tion of the tunnel spray cloud, can be used for a comparison between the 
experimental results obtained herein and the theoretical results . 

The modified inertia parameter is the independent variable 

chosen as the basis for comparing the experimental limits of impingement 
^max theoretical (fig. 15). For the Joukowski 0015 airfoil 

(fig. 15(a)) in the range of agreement of 

the experimental impingement limits on both upper and lower surfaces with 
those of theory is obtained for angles of attack of 0° and 4°. Fair 
agreement is obtained for Kq greater than 0.1. In this higher 

^,max the visual and colorimetric determinations of the experi- 

mental impingement limit, particularly on the lower surface, are more 
difficult than for Kn,Tnaw than 0.1 and may account for the lesser 

agreement with theoretical values when Kq is greater than 0.1. For 

the MCA 65]_-212 airfoil at a of 4° (fig. 15(b)) there is poor agree- 
ment on both upper and lower surfaces between experiment and theory, the 
theoretical limits being twice the experimental. Even on the upper sur- 
face, where at 4^ angle of attack the experimental limit is well defined, 
large discrepancies occur. As yet there is no reasonable explanation for 
these differences . 

The experimental local collection efficiencies ^ as a function of 
s are in good agreement with those obtained from theory for the Joukowski 
0015 airfoil__at both 0° and 4° angles of attack (fig. 16). A similar com- 
parison of |3 attempted for the M.CA 653_-212 airfoil at 4° angle of 

attack yielded very erratic results. The erratic nature of these data 
may be due in large part to the inconsistencies in the theoretical tra- 
jectory data. For a given impingement condition (Roq and K), total 

collection efficiency obtained by the tangent trajectory data of 
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reference 2 differs by more than 30 percent, for example, from deter- 

mined from an integration of the local collection efficiency results of 
the same reference (see eq. (4) herein). 

_ Theoretical and experimental values of total collection efficiency 
for the Joukowski 0015 at angles of attack of 0° and 4° and the NACA 

65]_-212 at an angle of attack of 4° are compared in figure 17 . For the 

Joukowski 0015 (figs. 17(a) and (b)) good agreement (±10^) between theory 
and experiment is obtained over the entire Kq range studied. The 

Ejjj value computed theoretically for the NACA 653_-212 airfoil (fig. 17(c)) 

is 25 and 100 percent greater than that obtained experimentally for 
^0,med v^a-lues of 0.09 and 0.008, respectively, because of the aforemen- 
tioned discrepancies in s^g^ and p. Similar comparisons between ex- 
periment and theory for the NACA Sbg-Olb airfoil at a of 4° (theoreti- 
cal in ref. 5) show the theoretical values to be 10 to 20 percent higher 
than the experimental values in the Kq range tested. 

_The___agreement between the theoretical and experimental evaluations 
^max considered good for the Joukowski 0015 and sat- 
isfactory for the NACA 652“015. No such agreement was obtained for the 

NACA 65^-212 airfoil, as previously discussed. These three airfoils are 

the only ones available for comparison of experimental and theoretical 
impingement values at this time. 


Correlation of Effect of Airfoil Sweepback 

The experimental impingement data substantiate the method of refer- 
ence 14 for predicting the impingement on a swept airfoil from data ob- 
tained on an unswept airfoil where (l) the wing can be considered two- 
dimensional or has a high aspect ratio and (2) the airfoil section in a 
plane perpendicular to the leading edge of the swept airfoil is the same 
section as that of the unswept airfoil. The application of the method 
of reference 14 to .the experimental impingement data presented herein is 
discussed in_appendix C, Typical experimental values of local collection 
efficiency 3 as a function of s for the 35° swept (NACA 653_-212 and 

6 ^ 1 ” 206 ) airfoils are presented for angles of attack (referenced to free- 
stream velocity direction) of 0° and 4.3° in figures 18(a)_and (b), re- 
spectively. The faired lines of figure 18 represent the 3 values cal- 
culated from the unswept experimental data by the method of reference 14. 
Good agreement between calculated swept and experimental swept values of 
3 was obtained for the angles of attack and the airfoils studied. Sim- 
ilar good agreement was obtained for values of and eL.. 

hicla. ni 


3689 



MCA TM 3839 


17 


Application of Tunnel Impingement Data to Flight 

Correlation of tunnel cloud characteristics -with those report ed in 
the literature for natural icing clouds . - Most of the reported data con- 
cernihg droplet size and liquid-water content of natural icing clouds 
(refs. 18 to 20) have "been obtained with rotating multicylinders that 
were permitted to ice for a known time interval. The rate of ice collec- 
tion on various size cylinders is matched with the theoretical collection 
of these cylinders in a manner that determines the droplet size and 
liquid-water content of the cloud (ref. 1). Similarly, the dye catch on 
various size cylinders can be matched to the theoretical cylinder catch 
to evaluate the tunnel dyed-water spray cloud properties. Another method 
of using the theoretical cylinder data to determine cloud properties is 
that of reference 15, wherein the droplet sizes are determined from dye- 
tracer impingement rates on a single stationary cylinder or body for 
which theoretical trajectories are available. A modification of the 
method in reference 15 is the use of a 36.5-percent symmetrical Joukowski 
airfoil. Details of the dye-tracer droplet-size analysis using the 36.5- 
percent Joukowski are given in appendix B. The liquid-water content of 
the dyed-water spray cloud was obtained by an aspirating tube (ref. 15). 

For droplets with diameters greater than 12 microns, the multi - 
cylinder matching, the single -cylinder solution, or the Joukowski airfoil 
solution each yield nearly the same absolute values of droplet size. 

Total liquid-water content, as measured by the aspirator, is nearly the 
same as that indicated by the multicylinder matching technique, and thus 
the VO Ixime -median droplet size is substantially the same for both the 
Joukowski -aspirator or multicylinder matching methods. The tunnel cloud 
properties and impingement data reported herein are based on the 
Joukowski-aspirator method. The difference or relation between the multi- 
cylinder matching and the Joukowski-aspirator evaluation of the tunnel 
spray cloud is illustrated in figure 19, which is a cross plot of the 
cloud properties as calibrated by the two techniques. The development of 
figure 19, a discussion of the relatively minor differences obtained, and 
the reasons for preferring the Joukowski-aspirator results over the multi- 
cylinder matching results are discussed in appendix B. 

To apply the experimental data herein to flight conditions, a pro- 
cedure is suggested and illustrated by a hypothetical problem, the condi- 
tions of which are as follows; 

(1) Meteorological design conditions, based on multicylinder data; 

Cloud volume-median droplet diameter, 15 microns 
Cloud total liquid-water content, 0.5 g/cu m 
Cloud droplet-size distribution, Langmuir ”D" 
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(2) Section characteristics: 

Airfoil section, MCA 653_-212 

Airfoil chord length, 10 ft 
Airfoil angle of attack, 4° 

(3) Operating conditions: 

True airspeed, 300 mph (261 knots) 

Pressure altitude, 10,000 ft 
Static air temperature, 12° F 

For the example, it is desired to determine the local water collec- 
tion rate at s^ of 0.02 and the lower surface limit of impingement 

^2, max’ 

The meteorological conditions are converted from the multi cylinder 
values to the Joukowski-aspirator values by using figure 19. The conver-,, 
sion results in a volume-median droplet size of 14.8 microns (Joukowski- , 
aspirator value) and a water-content ratio of the aspirator to the multi- 
cylinder match of 1.12. This ratio yields an aspirator total liquid-water 
content w^ of 0.56 g/cu m. Values of K and Roq are then calculated 

(0.03242 and 111.3, respectively). With this Rsq, a X/kg ratio of 0.325 
is obtained (fig. 6). The Kq calculated from K and X/Xg then 

amounts to 0.01054. In order to obtain ^ at the desired lower-surface 
location, s’ and s" must be obtained. The value of s" obtained from 
figure 10(f) and (g) at a of 0.01054 is 0.0037. From equation 

(5), s’ is then calculated to be 0.0163. The value of P is now deter- 
mined from figure 9(f) and (g) using Wie curves for 4® angle of attack and 
^0,med 0.01054; the result is a p of 0.14. The local collection rate 

at s^ of 0.02 is calculated from equation (2), which gives a value of Wp 
of 7.74 pounds per hour per square foot. 


In order to determine the limit of impingement on the lower surface 
of the airfoil, the maximum droplet size in the cloud droplet distribution 
must be established. For the Langmuir "D" distribution as well as the 
tunnel spray distribution, an average ratio of ^ax/^ed exists. 

Hence, the maximum droplet size for the example is 48 microns. A value 
of is now determined in a manner similar to that used to obtain 

^0,med values of ^/^s,max determined). 

The result of these calculations is a value of 0.0692, from 

which a value of s^ of 0.12 is obtained by use of figures 11(f) 

and (g). 
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Effect of droplet-size distribution on impingement characteristics « - 
A natural icing cloud may contain a distribution of droplet sizes (perhaps 
distribution types "B" to ”E" defined in ref, 12), or the droplets in 
a cloud may be all the same size ("A” distribution) . Many current design 
specifications are based on this assumption of an ”A" distribution. For 
large bodies, such as root sections of tapered wings or radomes (low 
Ko, range), the wider droplet-size distributions will result in higher 

values of Ej^ than will an "A” distribution. Therefore, icing -protection 

equipment designed for large bodies and an "A” droplet-size distribution 
may underestimate E^j and prove inadequate for some icing conditions. 

On the other hand, in the high Kq range typical for tip sections of 

tapered wings, helicopter blades, and instruments, the assumption of an 
"A” distribution may overestimate Ejj^ slightly when compared with an as- 
sumption of a wider droplet-size distribution. 

Limit of impingement is a function only of the maximum droplet size 
present in the distribution. Because typical distributions often contain 
droplets 2 to 3 times larger than the uniform size of the ”A" distribu- 
tion, the extent or limit of impingement will be markedly increased if a 
droplet-size distribution other than "A” is experienced. In addition, 
the 3 profile will be altered by different droplet-size distributions. 

Therefore, a droplet -size distribution that occurs relatively fre- 
quently in nature should be considered in the design of all icing- 
protection equipment. According to references 18 to 20, typical size 
distributions in nature range from a Langmuir "C” to "E". As previously 
discussed, the tunnel distribution of droplet sizes approximates a 
Langmuir "D”. Consequently, the droplet-size distribution inherent in 
the data reported herein is typical of that in many natural icing clouds, 
making these data suitable for design purposes. 


CONCLUDING REMARKS 

In an effort to obtain a general solution to the impingement char- 
acteristics of various airfoils, some investigators have suggested that 
correlation of impingement characteristics could be obtained if the body 
dimension used in the independent impingement parameters were based on 
airfoil thickness (ref. 21) or projected frontal height of the airfoil 
(ref. 22). In both of these references relatively good correlation was 
obtained for a limited number of airfoils at an angle of attack of 4°. 
Subsequent data obtained at other angles of attack, particularly 0*^, tend 
to show that the correlation at 4° was fortuitous. Much of the available 
theoretical E„^ data at 0° angle of attack is shown in figure 20(a) as a 

function of Kq^^,, where Kq,F based on the projected height of the 



20 


NACA TN 3839 


airfoil rather than the chord. (At 0° angle of attack the projected 
height is equal to the airfoil thickness.) It is apparent that at values 
of greater than 0.2 a wide deviation in total collection efficiency 

occurs, and the thick airfoils (thickness ratios 15 to 36,5^) have a 
higher collection efficiency than the thin airfoils (thickness ratios 4 
to 6^). Conversely, in the low ranges of the opposite trend is 

noted. 

A similar plot of data at a 4° angle of attack in figure 20(h) shows 
the same good correlation of data for the various airfoils as noted in 
reference 22. The apparent good agreement at the 4*^ angle of attack may 
he due to the fact that at a of 4° compared with a of 0° the projected 
frontal height of the thin airfoils more nearly approaches that of the 
thick airfoils (see fig. 5), thereby reducing the effect of the thickness 
ratio on impingement characteristics and parameters. The experimental 
data for collection efficiency show trends similar to figure 20 when 
plotted in the Kq y form. 

At present, there is no known parameter that accurately correlates 
all the available two-dimensional airfoil impingement data over realistic 
ranges of the independent variables. 


SUMMARY OF RESULTS 

The impingement characteristics of several airfoils obtained experi- 
mentally using a dye-tracer technique yield the following results: 

1. In general, the data show that the local and total water catch 
and the limit of impingement of airfoils are primary functions of the 
modified inertia parameter (in which airspeed and droplet and body size 
are the most significant variables) and airfoil thickness ratio. In 
addition, the local water collection rate and the extent of impingement 
on the airfoil surfaces depend on the airfoil angle of attack. Secondary 
factors affecting airfoil impingement characteristics are airfoil shape 
(for a given thickness ratio), small camber, and sweep angle. 

2. With an increase in the modified inertia parameter, the total and 
local collection efficiencies and the impingement limits also increase. 
For those airfoils of a comparable series operating at a typical flight 
value of the modified inertia parameter, a thickness ratio of 6 percent 
had total collection efficiencies of 1.7 and 2 times those of a 12- and 
16 -percent -thick airfoil, respectively. Airfoils with relatively blunt 
leading edges (Joukowski 0015) had higher total collection efficiencies 
than those with sharp leading edges (low-drag airfoils such as the NACA 
652-015), although the impingement limits for the sharper airfoils were 

greater than those for the blunt airfoils. 
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3. The experimentally determined local and total collection effi- 
ciencies and impingement limits for the Joukowski 0015 and KACA GSg-OlS 
airfoils are in good agreement with the theoretical values. No such 
agreement is obtained for the NACA 65]_-212 airfoil, 

4. Over the range of conditions studied, the experimental data sub- 
stantiate a previous method of predicting the impingement characteristics 
of swept airfoils (design section laid out perpendicular to the leading 
edge) from data for the unswept design section. 

5. Because of the typical droplet-size distribution of the tunnel 
spray, and the correlation of data by means of the modified inertia pa- 
rameter, the experimental results herein may be applied over a wide range 
of flight conditions. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 13, 1956 
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APPENDIX A 


SYMBOLS 


area,, sg ft 

volume of distilled water used to dissolve dye from blotter 
segments, ml 

percent concentration by weight of dye in water solution used in 


spray 


system, ^ 100% = x 100^ 

lb solution T.ro+o-n 


lb water 


airfoil chord length, ft 
cylinder diameter, in. 

droplet diameter, microns (3.28x10"^ ft) 

total collection efficiency in clouds of uniform droplet size, 
defined by eq. (3), dimensionless 

frontal height of airfoil projected parallel to free-stream veloc- 
ity direction divided by chord length, dimensionless 


inertia parameter, 


8.77xlO"^^P(id%Q 


, dimensionless 


Kq modified inertia parameter, ™ K, dimensionless 

P concentration of solution obtained from blotter segments, mg 
dye/ml solution 

RCq free-stream Reynolds number with respect to droplet, 

4.81xiO“^dpUo , 

= 'y'Kcp, dimensionless 

5 distance along surface referenced from zero-chord point divided 
by chord length, dimensionless 


exposure time, sec 

free-stream velocity, mi/hr or knots x 1,15 

local velocity at outer edge of boundary layer, mi/hr or 
knots X 1.15 
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total vater impingement rate in cloud of uniform droplet size, 
lb/(hr)(ft span) 

Wo local water impingement rate in cloud of uniform droplet size, 
l'b/(hr) (sq ft) 

w cumulative liquid-water content contained in droplets of sizes 
from d^ngy to any particular droplet size, g/cu m 

w-t total liquid-water content of cloud, g/cu m 

X distance along chord line from zero-chord point divided hy chord 

length, dimensionless 

y distance perpendicular to chord line divided by chord length, 
dimensionless 

a airfoil angle of attack, deg 

p local collection efficiency in cloud of uniform droplet size, 
defined by eq. (2), dimensionless 

X sweep angle, deg 

6 cylinder central angle, deg 

X true range of droplet as projectile injected into still air, ft 

Xg range of droplet as a projectile following Stokes' law, ft 

[1 viscosity of air, lb/(ft)(sec) 

p density of air, Ib/cu ft 

density of droplet, 62.4 Ib/cu ft 

0.423p^UqC 

(p independent impingement parameter, , dimensionless 

Subscripts: 

F frontal, projected parallel to free-Stream-velocity direction 

I lower surface 


max 


maximum 
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med volume -median 

n normal plane 

s surface 

u upper surface 

X coordinate parallel to free-stream-velocity direction 
Superscripts : 

“ weighted value due to effects of more than one droplet size 
’ referenced to surface location of Pxnax 

" location of P^ax 
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Air-water 

Volume-median 

Liquid-water 

pressure 

droplet diameter. 

content. 

ratio 

^ed^ 

^t^ 


microns 

g/cu m 

0.5 

21.2 

0.54 

.6 

16.8 

.45 

.8 

11.6 

.25 


Single Cylinders or 36.5-Percent Joukowski Airfoils 

In reference 15 a method for determining droplet-size distribution 
is described in which the theoretical trajectory results for cylinders 
(ref. l) are applied to experimental impingement obtained on cylinders 
with the dye-tracer technique. Reference 15 shows that the experimental 
pressure distribution about cylinders deviates considerably from theory; 
these surface pressure differences probably reflect ’unknown differences 
in the flow field ahead of cylinders and hence droplet trajectories about 
cylinders. In order to eliminate or reduce these unknown effects, a 36.5- 
percent-thick symmetrical Joukowski airfoil (coordinates listed in table 
II) has been selected herein for which the experimental pressure and thus 
velocity distributions are in good agreement with the theoretical values 
(shown in fig. 22) . In addition, a bluff configuration like the 36.5- 
percent Joukowski lends itself to accurate determination of the point 
where a droplet impinges on the surface, a critical factor in evaluating 
experimentally as well as theoretically the pertinent impingement variables 
®max^ {3, and E^. The 36.5-percent Joukowski was studied theoretically 

(unpublished trajectory data) -and with the dye- tracer technique. For the 
36.5-percent Joukowski experimental studies, airfoils of 5.47- and 16.32- 
inch chord were used. These chord sizes gave about the same leading-edge 
diameters as the 2- and 6- inch cylinders of reference 15 and thereby pro- 
vided a dimensional similarity for comparing the results from the two 
types of bodies . 

The experimental techniques and methods of determining droplet size 
from 36.5-percent Joukowski dye traces are identical to those detailed 
for cylinders in reference 15. In the Joukowski analysis, the surface 
distance from the zero-chord point is denoted as s, whereas in the 
cylinder analysis surface distance is given by the central angle 9 . 

Curves of s ^^^ as a function of K and <P , and (B as a function of 

Smax Sind s are given for the 36.5-percent Joukowski airfoil in figures 

23 and 24, respectively. These theoretical Joukowski results were ob- 
tained at the Lewis laboratory with the mechanical analog described in 
reference 23. Figures 23 and 24 for the 36.5-percent Joukowski airfoil 
are comparable to figures 15 and 16, respectively, of reference 15 and 
are used with the experimental results (unpublished) in the same manner. 
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APPEKDIX B 


DETERMIMTION OF CLOUD DROPLET SIZE AM) LIQUID-WATER 
CONTENT FROM MEASURED IMPINGEMENT RATES ON BODIES 

Cloud droplet size and liquid-water content can "be obtained from 
measured impingement rates on a body for which theoretical trajectory 
data are known (ref. 15). Most of the published cloud characteristics 
have been obtained by collecting ice on various size cylinders and match- 
ing these data to theoretical collection rates (refs. 18 to 20) . Similar- 
ly, the dye catch on various size cylinders can be matched with theory to 
evaluate the tunnel dyed-water spray cloud. A modification to the multi- 
cylinder technique is that of reference 15, wherein the impingement rates 
obtained by dye traces on one stationary cylinder suffice. A refinement 
to the method of reference 15 is the use of an airfoil (36.5^ symmetrical 
Joukowski herein) instead of a cylinder. A discussion and comparison of 
these methods follow. 


Multi cylinder Matching 

To calibrate the tunnel dyed-water spray cloud by the multicylinder 
matching technique, seven nonrotating cylinders (diam. of l/8, l/4, l/2, 

1, 2, 4, and 6 in.) were covered with absorbent material and separately 
exposed in the tunnel cloud. The total water catch (as measured by the 
dye collected) per unit time, frontal area, and velocity (= E^w^) is 

plotted in figure 21 as a function of cylinder diameter including thickness 
of absorbent material. The log- log plot of figure 21 is the conventional 
presentation for analyzing multicylinder data (ref. l) . By the matching 
method described in reference 1, the tunnel data of figure 21 can be 
matched to a theoretical Langmuir "D” droplet-size distribution. As dis- 
cussed in reference 1, there is usually considerable latitude in selecting 
the best theoretical fit to any multicylinder data. A Langmuir "D" dis- 
tribution is selected for the tunnel data because it provides for the air- 
foils studied herein, the best over-aJ.1 agreement of limits and rates of 
impingement with those calculated theoretically. This agreement was pre- 
viously discussed and illustrated in figure 16, where experimental and 
theoretical impingement rates (p) on the Joukowski 0015 airfoil are com- 
pared. Matching the tvinnel multicylinder data to a Langmuir "D’* droplet- 
size distribution yields the following spray cloud characteristics: 
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From such, a procedure dimensional plots of droplet diameter as a function 
of cximulative liquid-water content can he obtained. 


Comparison of Multicylinder, Single -Cylinder, 
and 36.5-Percent Joukowski Airfoil Solutions 

A curve of droplet diameter against cumulative liquid-water content 
based on data obtained from the Joukowski airfoils for each of the three 
tunnel spray conditions used herein is shown in figure 25(a) . In ad- 
dition to the 36.5-percent Joukowski solution- of droplet sizes, the 
multicylinder matching the Langmuir "D" solution obtained from figure 21 
and the single -cylinder solutions of reference 15 are presented for com- 
parison. The average solution from each of the three methods is repre- 
sented by the faired curves of figure 25. These averages are adjusted 
to common liquid-water contents. This adjustment is necessary for com- 
parison on a dimensional basis, because data for the three methods of 
resolution were not obtained from the identical array (number and spacing) 
of spray nozzles. It is apparent from the curves that the solutions are 
in reasonable agreement. 

Further comparisons of the methods of determining droplet size are 
presented in figure 25(b). For convenience, the cylinder size results 
of reference 15 are reproduced together with the 36.5-percent Joukowski 
results for the two chord sizes studied. The 36.5-percent Joukowski 
airfoil solutions in figure 25(b) are considered first, and each chord 
size gives a slightly different droplet-size distribution. However, the 
consistent body- size trend for droplet diameters less than 16 microns 
noted in the cylinder data (ref. 15) no longer exists. In addition, the 
over-all spread in droplet size at a particular liquid-water content value 
w is markedly reduced for droplet diameters greater than 16 microns . 
Ideally, different body sizes exposed to the same cloud should indicate 
the same droplet-size distribution. However, as discussed in reference 
15, the consistent body-size trend for cylinders may be the result of 
nonideal flow about cylinders or unaccounted-for droplet drag increments 
due to droplet acceleration. Furthermore, difficulties encountered in 
calculating accurately the theoretical trajectories in the range of low 
K values (less than about 0.7) also may contribute to the cylinder size 
trend when the theoretical data are applied in the analysis of the ex- 
perimental data. The air flow about the 36.5-percent Joukowski is nearly 
ideal, as previously discussed and illustrated in figure 22. This im- 
proved air flow and an increased accuracy of the theoretical trajectory 
data in the low K range for the 36.5-percent Joukowski airfoil evidently 
yield results nearer the ideal than those from cylinders, as shown in 
figure 25(b) . 
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In s-ummary, the 56.5-percent Joukowski solution of droplet size is 
preferred over a single or multicylinder solution "because (l) the con- 
sistent "body-size trend noted for cylinders is absent, and (2) the "body- 
size spread throughout the droplet-size range is reduced. An aspirating 
tube (ref. 15) is the preferred instrument in measuring total liquid- 
water content, because its collection efficiency is adjusted to be 100 
percent whereas the calculated efficiency for cylinders may be inaccurate 
in the low K range. 

Joukowski- aspirator values of volume-median droplet size and liquid- 
water content comparable to those previously tabulated in this appendix 
for the multicylinder matching method are as follows: 


Air- water 

pressure 

ratio 

Volume-median 

droplet 

diameter, 

"^ed^ 

microns 

Liquid-water 

content, 

^t^ 

g/cu m 

0.5 

18.6 

0.60 

.6 

16.7 

.47 

.8 

11.5 

.30 


A dimensionless droplet-size distribution as obtained by the Joukowski- 
aspirator method is discussed in the text and is presented in figure 3. 

Although the Joukowski -aspirator method of determining droplet size 
and liquid-water content of the dyed spray cloud is used to present the 
airfoil impingement characteristics herein, most of the published data 
on cloud characteristics have been obtained by the multicylinder matching 
method. The relation between the two methods for the range of conditions 
studied is obtained by a cross plot of the droplet sizes and liquid-water 
contents obtained by the two techniques. This cross plot, consisting of 
the data tabulated in this appendix, is shown as figure 19, and its use 
discussed in the text. 


"Values of Total Liquid-Water Content Reported in Reference 15 

The aspirator values of total liquid-water content as reported in 
reference 15 are in error because of an undetected recirculation of the 
dyed spray cloud. The error was incurred by operating the aspirator 
for a longer time than required for the air and dyed droplets in the tun- 
nel to recirculate. This phenomenon resulted in aspirator values of 
total liquid-water content higher than the true value by the amo\int of 
recirculated dye. Recent studies have evaluated the effect of recircu- 
lation for all spray conditions, and the effect is only" significant for 
values of total liquid- water content (and thus by definition volume- 
median droplet size) . The corrections to liquid- water content and 
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voltune -median droplet size (the latter based on analysis of cylinder im- 
pingement data) given in reference 15 are tabulated as follows : 


Air- water 

pressure 

ratio 

Liquid-water 
content , 

g/cu m 

Volume-median 
droplet diameter, 

^ed^ 

microns 

Refer- 

ence 

15 

Refer- 
ence 15 
cor- 
rected 

Refer- 

ence 

15 

Refer- 
ence 15 
cor- 
rected 

0.5 

0.70 

0.46 

14.8 

20.4 

.6 

.58 

.37 

12.0 

16.6 

.8 

.43 

.22 

7.6 

12.0 


The cylinder impingement, techniques, and method of solution for 
droplet size are unchanged from those reported in reference 15. The 
cylinders were exposed to the dyed spray cloud and then withdrawn from 
the tunnel before the circuit time of the tunnel air was complete (approx. 
14 sec at an airspeed in the test section of 152 knots) . The recirculated 
spray cloud for all conditions studied contains droplets less than 5 mi- 
crons in diameter with recirculation dying out completely in a few minutes. 
Even for the 13- inch-chord airfoils, these small droplets impinge in a 
very narrow band (less than l/l6 in.) at or near the leading edge. The 
effect of these recirculated droplets on all the airfoil impingement 
characteristics is considered negligible. Total liquid-water contents 
for the airfoils studied herein are measured with aspirator exposure 
times significantly less than the tunnel air circuit time. 
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APPENDIX C 


PREDICTION OF IMPINGEMENT ON SWEPT AIRFOIL 

FROM RESULTS OBTAINED ON UNSWEPT AIRFOIL 

Reference 14 presents a method for predicting the impingement on a 
swept wing from impingement data for an unswept airfoil section if the 
nnswept airfoil section is in the plane perpendicular to the leading edge 
of the swept wing. Limit of impingement and local and total collection 
efficiency are first determined with respect to flow conditions in the 
normal plane of the swept wing and then by geometry into the free-stream 
plane - 

The 35° swept NACA 65i-212 and 65i-206 models studied herein were 
layed out with the sections in a plane perpendicular to the leading edge. 
Data from these swept airfoils are therefore directly comparable with 
the data from the unswept airfoils modified by the method of reference 
14. Application of the unswept data as discussed in reference 14 to the 
swept airfoils is as follows (the NACA 65]_-212 airfoil is used as an 

example) : 

(1) For a particular value of Kq 3 is obtained from figure 

9(f) and (g) for several values of s* and angles of attack a. By 
using figure 10 and equation (5)_, s' is converted to s. The Kq 

value used to enter figure 9(f) and (g) is smaller than that at which the 
data were obtained. The value of Kq is smaller because Uq ^ is 

replaced by Uq^^^, where Uq^j^ = Uq^^i y; this has a greater effect on 
^,med accompanying increase in X/Xg (a function of Rcq and 

also decreased by cos r). . 

(2) The values of 3 obtained from step (l) are plotted against 
angle of attack a for constant values of s. 

(3) The p values are read from the plot__^described in step (2) at 
Ojj = a/cos y for several values of s. The p values thus obtained 

((Bj^) are in terms of a plane perpendicular to the leading edge of the 
swept airfoil. 

( 4 ) The Pn values from step (3) are multiplied by cos y to obtain 
values on the swept airfoil. Therefore, is in terms of a plane 

perpendicular to free-stream- velocity direction. . This latter definition 

of 3 (actually p^) is identical to the conventional one for unswept 
bodies . 
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( 5 ) The s (actually values of steps (l) to ( 3 ) are in terms of 

a plane perpendicular to the leading edge of the swept airfoil. These 
Sj^ values axe converted to corresponding s^^ values from the geometry 

of the "stretched" airfoil (the NACA 65]_-212 section in plane perpendicular 

to leading edge stretched to a thinner section in the free-stream plane) , 
For the MCA 65]_-212 and 652_-206, the difference between s^ and s^^. is 

of secondary importance. 

Unswept experimental ^ values modified by the preceding procedure 
are compared with experimental swept data in figure 13 and discussed in 
the text. 
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TABLE I. - IMPINGEMENT CHARACTERISTICS OF AIRFOILS 


Air- 

s 

Joukowskl 0015 


NACA 63p-015; 


water 


chord, 13 in. 


chord, 13 

in. 


pres- 



















sure 


Local efficiency 

P, for angle 

of attack, a. 

of - 

ratio 











OO 

2° 

4° 

8° 

0° 

4° 

8° 

12° 


0.5 

-0.08 

0.070 

0.035 

0.014 

<0.01 

0.044 

<0.01 




-.06 

.133 

.087 

.040 

<.01 

.089 

.014 




-.04 

.265 

.216 

.140 

.026 

.170 

.056 




-.02 

.540 

.547 

.477 

.235 

.371 

.212 

0.035 

0.000 


-.01 

.657 

.663 

.603 

.392 

.518 

.432 

.212 

.042 


0 

.688 

.733 

.687 

.547 

.770 

.735 

.525 

.335 


.01 

.657 

.697 

.673 

.605 

.518 

.629 

.581 

.546 


.02 

.540 

.572 

.575 

.590 

.371 

.511 

.546 

.588 


.04 

.265 

.350 

.377 

.461 

.170 

.335 

.410 

.476 


.06 

.133 

.217 

.253 

.353 

.089 

.210 

.294 

.364 


.08 

.070 

.127 

.165 

.262 

.044 

.140 

.224 

.301 


.10 

.034 

.072 

.104 

.195 

.021 

.115 

.182 

.252 


.15 

.20 

.010 

<.01 

.021 

<.01 

.037 

.017 

.081 

.047 


.056 

.027 

.115 

.056 

.168 

.115 




m 

0.378 

0.400 

0.392 

0.350 

0.302 

0.354 

0.340 

0.315 


su , max 

-0.189 

-0.167 

-0.117 

-0.075 

-0.199 

-0.106 

-0.029 

-0.032 . 


® 1 , max 

.189 

.240 

.292 

.435 

.199 

.350 

.453 

.623 

0.6 

-0.08 

0.045 


<0.01 

<0.01 

0.030 





-.06 

.105 


.020 

<.01 

.066 

<0.01 




-.04 

.235 


.075 

<.01 

.144 

.015 




-.02 

-485 


.385 

.185 

.336 

.121 

0.021 



-.01 

.601 


.520 

.354 

.489 

.292 

.156 

0.029 


0 

.637 


.623 

.520 

.725 

.664 

.489 

.285 


.01 

.601 


.614 

.590 

.489 

.605 

.547 

.497 


.02 

.485 


.523 

.570 

.336 

.489 

.518 

.526 


.04 

.255 


.347 

. 445 

.144 

.307 

.389 

.437 


.06 

.105 


.225 

.336 

.066 

.179 

.281 

.336 


.08 

.045 


.133 

.240 

.030 

.109 

.204 

.262 


.10 

.15 

.020 

.010 


.077 

.022 

.175 

.070 

.015 

.072 

.037 

.161 

.088 

.223 

.132 




.20 

<.01 

<.01 

.032 

.011 

.049 

.015 






0.310 


0.317 

0.313 

0.252 

0.279 

0.295 

0.254 


s 

u.max 

-0.158 


-0.105 

-0.058 

-0.186 

-0.073 

-0.027 

-0.028 


tjraax 

.158 


.246 

.379 

.186 

.315 

.442 

.508 

0.8 

-0.08 

<0.01 

<0.01 

<0.01 

<0.01 

1 

0.008 





-.06 

.037 

.011 

<.01 

<.01 

.022 





-.04 

.115 

.051 

.017 

<.01 

.071 





-.02 

.380 

.287 

.207 

.057 

.212 

0.049 




-.01 

.490 

.415 

.355 

.201 

.367 

.183 

0.033 

0.017 


0 

.553 

.521 

.467 

.365 

.679 

.501 

.267 

.200 


.01 

.490 

.505 

.484 

.453 

.367 

.516 

.449 

.383 


.02 

.380 

.398 

.415 

.448 

.212 

.384 

.467 

.434 


.04 

.115 

.203 

.257 

.353 

.071 

.203 

.333 

.347 


.06 

.037 

.090 

.135 

.241 

.022 

.100 

.200 

.242 


.08 

<.01 

.037 

.055 

.150 

.008 

.049 

.133 

.175 


.10 

<.01 

.014 

.027 

.087 



.033 

.100 

.133 


.15 

<.01 

<.01 

<.01 

.025 


.011 

.033 

.059 


.20 

<.01 

<.01 

<.01 

<.01 



.017 

.017 


^m 

0.219 

0.213 

0.197 

0.198 

0.157 

0.169 

0.193 

0.156 


®u,max 

-0.095 

-0.069 

-0.054 

-0.032 

-0.107 

-0.034 

-0.019 

-0.020 


S 1 

i-jinax 

.095 

.119 

.154 

.234 

.107 

.192 

.280 

,319 
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TABLE I. - Continued. IMPINOEMENT CHARACTERISTICS OP AIRFOILS 


Air- 

water 

pres- 

3 

NACA 652-015; 
chord, 13 In. 


NACA 651-212 
chord, 13 In 

* 

sure 










ratio 


Local efficiency. 

F, for 

angle of attack, a. 

Df - 



0° 

4° 

8° 

12° 

0° 

20 

4° 

8° 

0.5 

-0.08 

0.039 

<0.01 



0.035 

0.014 

<0.01 



-.06 

.072 

.011 



.070 

.040 

.014 



-.04 

.157 

.038 




.154 

.084 

.035 



-.02 

.345 

.168 

0.039 

0.014 

.336 

.224 

.113 

0.021 


-.01 

.542 

.378 

.280 

.085 

.537 

.490 

.301 

.168 


0 

.821 

.749 

.638 

.329 

.757 

.770 

.699 

.546 


.01 

.542 

.629 

.651 

.582 

.336 

.505 

.531 

.560 


.02 

.345 

.480 

.562 

.550 

.197 

.308 

.349 

.453 


.04 

.157 

.294 

.413 

.440 

.084 

.168 

.197 

.322 


.06 

.072 

.182 

.301 

.344 

.037 

.097 

.154 

.252 


.08 

.039 

.126 

.238 

.270 

.014 

.070 

.115 

.201 


.10 

.024 

.104 

.191 

.224 

<.01 

.049 

:084 

.168 


.15 

<.010 

.070 

.126 

.154 

<.01 

.028 

.057 

.105 


.20 


.033 

.085 

.115 

<.01 

.014 

.028 

.070 



0.289 

0.339 

0.384 

0.316 

0.274 

0.326 

0.327 

0.344 


®u,max 

-0.280 

-0.095 

-0.035 

-0.027 

-0.250 

-0.176 

-0.109 

-0.036 


® 1 , max 

.280 

.422 

.527 

.600 

.202 

.350 

.460 

.580 

0.6 

-0.08 

0.021 

- 



0.029 

<0.01 

<0.01 



-.06 

.045 

<0.01 



.058 

.022 

<.01 



-.04 

.113 

<.01 



.121 

.058 

.015 



-.02 

.287 

.088 

0.023 

<0.01 

.307 

.205 

.073 

0.014 


-.01 

.467 

.248 

.117 

.044 

.482 

.464 

.226 

.117 


0 

.759 

.608 

.518 

.258 

.715 

.737 

.562 

.512 


.01 

.467 

.606 

.574 

.508 

.365 

.473 

.502 

.554 


.02 

.287 

.467 

.489 

..497 

.161 

.292 

.336 

.454 


.04 

.113 

.258 

.347 

.394 

.058 

.153 

.190 

.322 


.06 

.045 

.150 

.248 

.307 

.029 

.088 

.132 

.242 


.08 

.021 

.099 

.183 

.233 

.011 

.044 

.103 

.190 


.10 

.012 

.072 

.146 

.177 

<.01 

.029 

.073 

.146 


.15 


.043 

.072 

.117 

<■01 

.011 

.037 

,088 


.20 


.02g 

.044 

.072 


<.01 

.015 



.051 



0.216 

0.270 

0.275 

0.243 

0.237 

0.268 

0.267 

0.302 


^u, max 

-0.200 

-0.069 

-0.032 

-0.024 

-0.199 

-0.129 

-0.085 

-0.029 


S 1 

.200 

.373 

.473 

.520 

.170 

.288 

.394 

.534 


1 , max 






' 



o 

CD 

-0.08 






<0.01 





-.06 

0.013 




.017 

<0.01 




-.04 

.058 




.051 

.017 




-.02 

.195 

0.040 

0.01 

— 

.201 

.092 

0.025 



-.01 

.359 

.175 

.049 

<0.01 

.425 

.300 

.141 

0.025 


0 

.724 

.515 

.283 

.149 

.699 

.674 

.567 

.416 


.01 

.359 

.549 

.471 

.379 

.316 

.376 

.500 

.532 


.02 

.195 

.372 

.400 

.400 

.116 

.191 

.327 

.425 


.04 

.058 

.172 

.266 

.300 

.033 

.081 

.159 

.259 


.06 

.013 

.092 

.167 

.216 

<.01 

.033 

.092 

.175 


.08 


.049 

.116 

.149 

<.01 

.017 

,059 

.122 


.10 


.025 

.083 

.116 


.011 

.033 

.083 


.15 


.013 

.033 

.049 


<.01 

.013 

.033 


.20 


<.01 

.01 

.025 



<.01 

.017 


^m 

0.147 

0.167 

0.167 

0.143 

0.186 

0.186 

0.200 

0.191 


^u,max 

-0 . 120 

-0.037 

-0.024 

-0.012 

-0.121 

-0.070 

-0.041 

-0.015 


^ l,max 

.120 

.214 

.285 

.346 

.081 

.152 

.248 

.302 
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TABLE I. - Continued. IMPINGEMENT CHARACTERISTICS OF AIRFOILS 


Air- 

water 

pressure 

ratio 

s 

Joukowskl 0015; 
chord, 96 in. 

NACA 65,2-216; 
chord, 96 in. 

Local efficiency, p, for angle 
of attack, a, of - 

0° 

2° 

40 

0° 

40 

0.5 

-0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.02 

.03 

.04 

.06 

.08 

0.057 

.139 

.206 

.274 

.334 

.368 

.334 

.274 

.206 

.139 

.057 

.015 

<.01 

.0 

0.031 

.083 

.137 

.229 

.305 

.365 

.400 

.336 

.275 

.229 

.118 

.046 

<.01 

<.01 

0.015 

.031 

.061 

.121 

.213 

.320 

.383 

.397 

.350 

.294 

.187 

.095 

.020 

<.01 

0.031 

.077 

.152 

.282 

.395 

.442 

.371 

.243 

.131 

.059 

.015 

<.01 

<.01 

.0 

0.00 

<.01 

.015 

.036 

.164 

.317 

.466 

.437 

.334 

.268 

.146 

.077 

.031 

.015 


m 

0.092 

0.096 

0.098 

0.075 

0.078 

®u,max 
® 1, max 

-0.08 

.08 

-0.06 

.10 

-0.05 

.13 

-0.09 

.08 

-0.03 

.13 

0.6 

-0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.02 

.03 

.04 

.06 

.08 

0.015 

.078 

.141 

.203 

.266 

.297 

.266 

.203 

.141 

.078 

.015 

<.01 

.0 

<0.01 

.031 

.063 

.141 

.221 

.282 

'.311 

.263 

.203 

.141 

.053 

.015 

<.01 

0.0 

<.01 

.015 

.051 

.123 

.210 

.277 

.292 

.251 

.203 

.109 

.040 

<.01 

<.01 

<0.01 

.034 

.069 

.234 

.350 

.397 

.314 

.191 

.051 

.022 

<.010 

.0 

0.0 

<.01 

.020 

.058 

.266 

.360 

.354 

.275 

.203 

.088 

.031 

<.01 

<.01 


Em 

0.060 

0.063 

0.059 

0.055 

0.057 

s 

u,max 

Si 

1 ,max 

-0.055 

.055 

-0.04 

.07 

-0.03 

.09 

-0.07 

.04 

-0.02 

.10 

0.8 

-0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.02 

.03 

.04 

.06 

.08 

0.0 

.011 

.030 

.079 

.146 

.176 

.146 

.079 

.030 

.011 

.0 

0.0 

<.01 

.027 

.091 

.145 

.174 

.145 

.079 

.039 

<.01 

.0 

0.023 

.068 

.137 

.160 

.136 

.079 

.020 

<.01 

0.0 

.206 

.256 

.175 

.061 

<.01 

.0 

<0.01 

.066 

.198 

.221 

.152 

.089 

<.01 

<.01 


Em 

0.025 

0.025 

0.023 

0.026 

0.023 


®u,max 
® 1 , max 

-0.03 

.03 

-0.02 

.04 

<-0.01 

.05 

-0.01 

.02 

<-0.01 

.05 
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TABLE I. - Concluded. IMPINGEMENT CHARACTERISTICS OP AIRFOILS 


Alr- 

waterl 

pre3-| 

sure 

ratio! 


- 0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.020 

.03 

.04 

.06 

.08 




NACA 65^^-212; 
chord, 72 In. 


NACA 65j ^- 206 ; 
chord, 72 In. 


‘^NACA G5j ^- 212 ; 

chord, 87.9 In.; 
swept 35° 


^NACA 65i - 206 ; 
chord, 87.9 in. 
swept 35° 


0.038 

.085 

.187 

.311 

.492 

.599 

.382 

.170 

.058 

.024 

.013 

<.01 

.0 


2 ° 


< 0.01 

.020 

.062 

.230 

.438 

.592 

.507 

.322 

.164 

.098 

.038 

.014 

<.01 

<.01 


0.106 


Local collection efficiency, P, for a ngle of attack, a, 

2 ° 


of - 


O.OOO] 

<.01 

.019 

.0531 

.266 

.492 

.559 

.433 

.302 

.202 

.115 

.071 

.031 

.015 


0.122 


0 ° 


0.034 

.059 

.095 

.178 

.426 

.795 

.117 

.016 

<.01 

<.01 


0.158 


2 ° 


< 0.01 

<.01 

<•01 

.034 

.187 

.720 

.328 

.105 

.057 

.043 

.030 

.022 

.011 

<.01 


0.132 


0.0 

<.01 

<.01 

.034 

.620 

.486 

.276 

.189 

.147 

.098 

.067 

.037 

.027 


0.147 


0 ° 


0.017 

.040 

.092 

.215 

.365 

.472 

.290 

.129 

.042 

.019 

<.01 

<.01 

.0 


0.0 

<.01 

.0151 

.070 

.202 

.337 

.321 

.196 

.115 

.068 

.028 

.019 

<.01 

<.01 


0.077 


0.0 

<.01 

.015 

.096 

.270 

.415 

.337 

.242 

.170 

.098 

.064 

.030 

.018 


0.108 


0.013 

.024 

.040 

.099 

.294 

.605 

.075 

.021 

<.01 

<.01 

0 


0.127 


2 ° 


0.0 

<.01 

<.01 

.091 

.629 

.238 

.117 

.070 

.049 

.036 

.026 

.011 

<.01 


0.129 


< 0.01 

<.01 

.394 

.397 

.246 

.174 

.132 

.087 

.062 

.033 

.019 


0.128 


0.6 


u,rnax 

^l,raax 


- 0.10 

.07 


. 0.06 

.12 


- 0.03 

.18 


- 0.13 

.025 


- 0.04 

.16 


- 0.02 

.25 


- 0.08 

.06 


- 0.03 

.18 


- 0.02 

.13 


• 0.09 

.03 


- 0.02 

.20 


k-o.oi 

>.25 


- 0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.020 

.03 

.04 

.06 

.08 


0.011 

.043 

.115 

.262 

.426 

.527 

.312 

.108 

.031 

<.01 

.0 


< 0.01 

<.01 

.028 

.124 

.342 

.502 

.433 

.231 

.115 

.063 

.015 

<.01 

<.01 

<,01 


0.0 ' 
<.01 
.028 
.178 
.438 
.501 
.380 
.257 
.173 
.079 
.044 
.011 
<.01 


0.019 

.032 

.054 

.129 

.352 

.661 

.087 

<.01 

<.01 

,0 


0.00 

<.01 

<■01 

.019 

.152 

.671 

.280 

.080 

.037 

.027 

.018 

.010 

<.01 

<.01 


0.000 

.019 

.506 

.455 

.246 

.162 

.117 

.067 

.040 

.020 

.012 


|< 0.01 

.023 

.059 

.157 

.320 

.420 

.236 

.082 

.019 

<.01 

<.01 

.0 


0.0 

<.01 

.028 

.143 

.280 

.280 

.169 

.089 

.047 

.023 

.014 

<.01 

<.01 


0,0 

.037 

.174 

.336 

.277 

.196 

.123 

.066 

.038 

<.014 

<.01 


|< 0.01 

.015 

.024 

.070 

.249 

.615 

.052 

.013 

<.01 

.0 


0.0 

.059 

.555 

.221 

.087 

.049 

.037 

.027 

.019 

<.01 

<.01 


0.078 


0.077 


0.092 


0.118 


0.107 


0.il7 


0.072 


0.059 


0.068 


0.096 


0.100 


k 0.101 


u,max 

,raax 


- 0.06 

.03 


• 0.04 

.08 


- 0.02 

.13 


- 0.09 

.02 


- 0.03 

.12 


- 0,01 

.17 


- 0.05 
.04 


- 0.02 

.09 


- 0.01 

.14 


- 0.05 

.02 


- 0.01 

.14 


- 0.03 

-.02 

-.015 

-.010 

-.005 

0 

.005 

.010 

.015 

.020 

.03 

.04 

.06 

.08 


ko.oi 

<.01 

.038 

.146 

.328 

.433 

.193 

.025 

<.01 

.0 


0.00 

<.01 

.029 

.203 

.380 

.306 

.124 

.039 

<.01 

.0 


< 0.01 

<•01 

<.01 

.043 

.262 

.382 

.279 

.155 

.073 

.020 

<.01 

<.01 


< 0.010 
.016 
.025 
.091 
.324 
.645 
.022 
.01 


Em 


0.048 0.045 


0.050 


0.094 


0.000 

.063 

.624 

.241 

.042 

.018 

.014 

<.01 

<.01 


< 0.01 
<.01 
.383 
.415 
.197 
.106 
.065 
.038 
.023 
< .01 
.0 


0 

< 0.01 
< .01 
.071 
.215 
. 333 , 
.122 
,024 
<.01 
.0 


0.084 


0.077 


0.042 


0 

kO.Ol 
.01 
.071 
.221 
.237 
.118 
.043 
.016 
.01 
,01 


< 0.01 
.075 
.237 
.198 
.114 
.063 
.022 
<-01 
<.01 
< .01 


.01 

.01 

.016 

.030 

.193 

.545 

.028 

.0 


0.012 
.490 
.193 
.063 
.030 
.020 
.016 
.014 
< .010 


0.040 


0.038 


0.071 


< 0.01 
.324 
.383 
.232 
.155 
.111 
.068 
.045 
.019 
:.01 


- 0.01 

>.20 


0.217 

.312 

.166 

.091 

.054 

.029 

.016 

<.01 

<.01 


0.065 


u,max 
I ,max 


- 0.04 

.02 


- 0 . 02 
.03 


- 0.01 

.07 


- 0.05 

<.01 


- 0.01 

.05 


k-o.oi 
.08 


- 0.03 

.02 


- 0.02 

.05 


<- 0.01 

.09 


• 0.04 

.01 


k-o.oi 

.07 


<- 0.01 
.14 


^Design section in plane normal to leading edge. 
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TABLE II. - COOEDIMTES OF 36.5- 
PERCEIIT-THICK JOUKOWSKI AIRFOIL 


X 

y 

s 

0.0000 

0.0000 

0.000 

.0024 

.0248 

.025 

.0095 

.0492 

.050 

.021 

.0726 

.078 

.0372 

.0945 

.104 

.0575 

.1146 

.132 

.0816 

.1325 

.161 

.1094 

.1479 

.194 

.1404 

.1606 

.227 

.1741 

.1704 

.262 

.1884 

.1734 

.275 

.2029 

.1760 

.291 

.2104 

.1771 

.299 

.2400 

.1810 

.328 

.2803 

.1818 

.368 

.2884 

.1817 

.376 

.3713 

.1752 

.458 

.4563 

.1594 

.548 

.5408 

.1369 

.633 

.6229 

.1106 

.723 

.7008 

.0834 

.800 

.7729 

.0580 

.879 

.8376 

.0364 

.948 

.8934 

.0199 

1.003 

.9388 

.0088 

1.049 

.9724 

.0027 

1.082 

.9931 

.0003 

1.103 

1.0000 

.0000 

1.110 


coor 
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Figure 1, - Typical installflrfcioiL of airfoil (MCA 653_-212) witk blotter 
attacbed, in 6- by 9- foot test section of icin^ tunnel. 




! 




t 

(d) 




t 

(J) 


^Designated section Is In plane perpendicular to 
leading edge; leading edge Is swept 35°. 

Figure 2. - Airfoil sections studied with dye-tracer technique. Sections laid out 
In free-stream direction. 
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Figure 3, - Dimensionless droplet-size distribution of tunnel spray obtained on 36.5- 
percent Joukowski airfoil (see appendix B), 
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Figure 4, - Typical blotter records from airfoils after exposure to dye 
cloud wltb. punched segments removed for colorimetric analysis. 
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Flgtcre 5. - Effect of angle of attack on ratio of projected height to 
chord length of mrious airfoils. (Reference is free-stream-velocity 
direction. ) 


'4 6 10 20 40 60 100 200 400 600 1000 

Free-stream Reynolds number, Roq 

Figure 6. - Droplet range ratio as function of droplet Reynolds nvimber (data from ref. 12). 
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Upper surface Lower surface 

Distance along airfoil surface, s 


Figure 7. - Typical water impingement rates on NACA 652-212 airfoil. Airspeed, 152 knots; volume- 
median droplet diameter, 16.7 microns; angle of attack, 4°. 
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Local collection efficiency 


Airfoil 



(c) Shape effect. Angle of attack, 0°,* 13-inch chord. 

Figure 8. - Concluded. Effect of airfoil geometry and angle of attack on local collection 
efficiency. Airspeed, 152 knots; volume -median droplet diameter, 16.7 microns. 
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Itodlfled inertia petrameter, 

(d) MCA eSg-OlS airfoil. Chord, 13 Inchea. 

Figure 9. - Continued. Local collection efficiency of airfolla as function of modified 
inertia parameter and surface location. 






Modified inertia parameter, Kq 

(e) NACA 65,2-216 airfoil. Chord, 96 inches. 

Figure 9. - Continued. Local collection efficiency of airfoils 
as function of modified inertia parameter and surface location 













Modified inertia parameter, K 


0,med 

(h) WACA 65^-206 airfoil. Chord, 72 inches. 

Figure 9. - Continued. Local collection efficiency of airfoils as function 
of modified inertia parameter and surface location. 
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Modified inertia parameter, Kq 


(i) NACA 65^^-212 airfoil, swept 35° (data in free-stream direction). 
Chord, 87.9 inches. 

Figure 9. - Continued. Local collection efficiency of airfoils as function 
of modified inertia parameter and surface location. 


NACA TN 3839 






Local collection efficiency 



MCA TE 3839 




Surface location of 


MCA TU 3839 









Surface location of 


58 


MCA TN 3859 




(i) NACA 65X-212; 55° sweep. (j) NACA 65x206; 35° sweep. 


Figure 10, - Concluded. Surface location of maximum local collection efficiency as 
function of modified inertia parameter (max. and foremost points tabulated as s 
distance. 
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(a) and (b) Joukowski 0015. (c) NACA 632-015* Chord, 13 Inches. 



(d) NACA 652-015. Chord, 13 inches. (e) NACA 65,2-216. Chord, 

96 Inches . 


Figure 11. - Limits of impingement on airfoils as function of modified inertia 
parameter. 






Limit of impingement 



1 .1 .01 .1 
Modified inertia parameter, 


(i) NACA 65;j^-212, 35° sweep. (j) NACA 65i-206, 35° sweep. 

Chord, 87.9 inches. Chord, 87.9 inches. 

Figure 11. - Concluded. Limits of impingement on airfoils as function of 
modified inertia parameter. 
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Angle of attack (geometric ) , a, deg 


(d) Chord, 96 inches; 0.0104 

Pigtrre 13. - Effect of angle of attack on experimental total collection 
efficiency of six airfoil sections. Airspeed, 152 knots; vol-ume- 
median drop diameter, 16,7 microns. 
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.1 .2 .3 .4 .5 .6 

Distance along airfoil surface, s 


(a) Joukowski 0015 airfoil; 96-inch chord. 

Figure 14. - Local velocity dlstrihutlon on airfoil. Airspeed, 152 
knots; pressure, 28.1 inches of mercury; temperature 50° F. 
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Diatance along airfoil surface, s 

(b.) HACA 65i**212 airfoil; angle of attack, 4°, 

Figure 14. - Concluded. Local velocity distribution- on airfoil. 
Airspeed, 152 knots; pressure, 28.1 Inches of mercury; tenpera- 
txtre, 50° F, 
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(a) Joukowski 0015 airfoil; angle of 
attack, 0° and 4°« 


Figure 15, - Comparison of experimental 



(b) NACA 65-j^-212 airfoil; angle of 
attack, 4P . 
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(a) Angle of attack, 0°. 



(b) Angle of attack, 4°. 


Figure 16. - Comparslon of experimental and theoretical local collection efficiency on 
15-percent symmetrical Joukowskl airfoil. Angles of attack, 0° and 4°; chord length, 
13 and 96 Inches; volume-median droplet size, 16.7 microns; airspeed, 152 knots. 
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Theoretical (ref. 5) 


96-In, chord 


(a) Joukowski 0015; angle of attack, 0°, 


Theoretical (ref, 5) 


o5 0 
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(h) Joukowski 0015; angle of attack, 4°, 


Theoretical (ref. 2)-i 

I I » t I 1 1 I L 


7 2- In. chord 


Modified inertia parameter, Ko,med 
(c) WACA 65;j_-212; angle of attack, A*' 


Figure 17. - Comparison of experimental and 
theoretical total collection efficiencies. 





Local collection efficiency. 


Airfoil 


NACA 653_-212 


MCA TN 3839 


Experimental swept 
Calculated from unswept 
esqjerimental data 


1 

NACA 65^ 
1 

— 

-206 












,! 


- 



04 .-.OS G. .02 .1 

Upper surface Lower surface 

Distance along airfoil surface, s 

(a) Angle of attack, 0°. 

Figure 18. - Comparison of local collection ef- 
ficiency for 35° swept airfoils with values 
calculated from unswept- airfoil data (method 
of ref. 14). Speed, 152 knots; volume- median 
droplet size, 16.7 microns. (All values in 
terms of free- stream direction. ) 
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Airfoil 


NACA 65-, -212 


NACA 65-, -206 


Experimental swept 
Calculated from unswept 
experimental data 


.06 .08 
Tipper surface Lower surface 

Distance along airfoil surface, s 

(b) Angle of attack, 4.3°. 

Figure 18. - Concluded. Comparison of local collection efficiency for 
35° swept airfoils with values calculated from unswept-airfoil data 
(method of ref. 14). Speed, 152 knots j volume- median droplet size, 
16.7 microns. (All values in terms of free-stream direction.) 
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Figure 20. - Correlation of total collection efficiency with inertia parameter modified 
by projected height. 
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Water collected ~ / 

“(Time) (Frontal area) (Velocity) " ^ 

Figure 21, - Water collected "by various cylinder sizes. Free- 
stream velocity, 152 knots; air temperature, 50° F; air 
pressure, 28.1 inches of mercury. 
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Figure 23. - Limit of impingement on upper or lower surface of 36 . 5-percent-thlck symmetrical 
Joukowski airfoil. Angle of attack, 0°. 
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Figure 24. - Local collection efficiency of 36 . 5-percent-thlck symmetrical 
Joukowski airfoil. Angle of attack, 0°. 
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Flgvjre 25, - Droplet-size distributions from experimental impingement 
rates on cylinders and 36.5-percent Joukowski airfoils. Comparison 
of methods of solution. 
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Figure 25. - Concluded. Droplet-size distributions from experimental impingement 
rates on cylinders and 36.5-percent Joukowski airfoils. Comparison of methods 
of solution. 
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